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Abstract 

We study the Drell-Yan production cross section and structure func- 
tions in proton (deuteron)-nucleus collisions using the Color Glass Con- 
densate formalism. The nucleus is treated in the Color Glass Condensate 
framework which includes both higher twist effects due to the inclusion 
of multiple scatterings and leading twist pQCD shadowing due to the 
small X resummation, while the proton (or deuteron) is treated within the 
DGLAP improved parton model. In particular, the Drell-Yan structure 
functions are used in order to investigate the Lam- Tung relation at small 
X, which is known to be identically zero at leading twist up to Next-to- 
Leading order, and is thus a good playground for studying higher twist 
effects. In agreement with this, we find that violations of this relation are 
more important for low momentum and invariant mass of the Drell-Yan 
pair, and also in the region of rapidity that corresponds to smaller values 
of X in the nucleus. 
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1 Introduction 



The recent data from deuteron-gold collisions in the forward rapidity region at 
RHIC [1] have caused a great excitement in the high energy heavy ion physics 
community. The observed suppression of the liadron production transverse 
momentum spectrum in deuteron-gold collisions as compared to proton-proton 
collisions (normalized to the number of binary collisions) was qualitatively pre- 
dicted by the saturation physics and the Color Glass Condensate formalism [2-5] 
and is now understood quantitatively using saturation based models [6 10]. 

While the produced hadron spectra in deuteron-gold collisions in mid ra- 
pidity can be described in terms of classical (Glauber like) multiple scattering 
[11-21], one needs to go beyond a simple multiple scattering picture in order to 
understand the suppression of the forward rapidity spectra. In the Color Glass 
Condensate formalism, this suppression is due to the evolution of the target 
nucleus wave function with x (or equivalently, with rapidity) described by the 
JIMWLK evolution equations [22-28,2-4] (or their large Nc limit known as the 
Balitsky-Kovchegov (BK) equation [29,30]). 

To further probe the extend to which saturation is the dominant physics 
in the forward rapidity region at RHIC, and in order to clarify the role of the 
additional effects, such as parton recombination, cold matter energy loss, etc., 
which have recently been suggested [31-33] as the reason for the suppression 
pattern, it is essential to consider the electromagnetic probes of the Color Glass 
Condensate, such as photon and dilepton production [34-36] or photon-hadron 
correlations [37] . The Color Glass Condensate formalism would predict a similar 
suppression pattern for production of photons and dileptons while the parton 
recombination models should not since photons and dileptons do not interact 
strongly. Therefore, observation of (or lack thereof) a similar suppression in 
the production of photons (or dileptons) as in hadron production would firmly 
establish saturation as the underlying physics of the phenomena in the for- 
ward rapidity region at RHIC. In this paper we consider dilepton production 
in proton-nucleus collisions and provide quantitative predictions for Drell-Yan 
structure functions and angular correlations in the kinematic region appropri- 
ate for the PHENIX experiment at RHIC and CMS and ATLAS experiments 
at LHC [38]. 

2 Drell-Yan cross-section in the CGC formalism 
2.1 Generalities 

Electromagnetic observables (such as photons and dileptons) are cleaner probes 
of new phenomena than hadronic ones since they do not interact strongly and 
are not sensitive to the non-perturbative, and thus poorly understood, physics 
of hadronization. However they do suffer from the fact that their production 
rates is much smaller than the hadron production rate, due to smallness of the 
electromagnetic coupling constant, and that in case of photons, measuring them 
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precisely poses new experimental difSculties. 

When looking for new physics effects, it helps enormously when the leading 
order processes which constitute the background vanish. This is the case for 
the so called Lam- Tung relation [39,40] in dilcpton production in pcrturbative 
QCD. It is known that the Leading Order (LO) pQCD partonic cross sections 
satisfy the Lam- Tung relation. This is the analog of the Callan-Gross relation 
between the structure functions Fi and F2 in Deeply Inelastic Scattering (DIS) 
which reads F2 = 2xFi or alternatively, Fl = F2 — 2xFi = 0. This is an 
indication that the photon probe scatters off of a spin half object (a fermion) 
in the target and as such, helped establish quarks as the constituents of the 
target proton. The Callan-Gross relation is, however, not exact and receives 
corrections due to radiation of gluons such that Fl ~ ag xG{x, Q^). 

The Lam- Tung relation on the other hand does not receive any ag corrections 
(the first correction is of order a^) and therefore can be a very sensitive probe 
of new (beyond standard pQCD) physics, such as higher twists effects and/or 
a change from DGLAP to BFKL kinematics signified by the change of the 
anomalous dimension of the gluon distribution function. Here, we use the Color 
Glass Condensate formalism, which includes both the multiple scattering and 
BFKL anomalous dimension into effect, in order to investigate the possible 
beyond the standard pQCD physics in the Lam- Tung relation. 

2.2 The cross-section and hadronic tensor 

In [34] and [35], the expressions for photon and dilepton production cross sec- 
tions in proton (deuteron)-nucleus collisions were derived, using a hybrid de- 
scription in which the proton is described by ordinary parton distributions while 
the nucleus is described in terms of classical color sources and the color field they 
produce. In such a description, the leading process for the emission of a photon 
(real or virtual) is a quark coming from the proton, scattering off the classical 
color field of the nucleus, and radiating a photon. Two terms contribute to the 
production amplitude, depending on whether the photon is emitted before or 
after the scattering of the quark off the nucleus, as illustrated in figure 1. There 



Figure 1: The two diagrams contributing to the emission of a photon in pA 
collisions. The black dot denotes the Wilson line that resums all the multiple 
scatterings of the quark off the color field of the nucleus. 

could potentially be a third diagram, where the emission of the photon occurs 
inside the nucleus, so that the quark can interact with the nucleus before and 
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after the photon emission. However, such a contribution is suppressed by the 
inverse of the colhsion energy due to the Lorentz contraction of the nucleus. 
If we denote p the momentum of the incoming quark, q the momentum of the 

outgoing quark, and k the momentum of the emitted photon, the amplitude for 
photon emission in a given configuration of the target color field reads [34,35] : 

M^{p\q,k) = 2wS{p--q--k-) j (fx± e'^-i^+'^^-P-^^"'^ ^V{x±)-1 

In this formula, the nucleus is assumed to be moving in the positive z direction. 
We have neglected the mass of the quark, and Cq is the electrical charge of 
the quark in units of the electron charge e. V{x±) is a Wilson line defined as 
follows : 



V{x±_) = T_exp 



-ig / dz A+{z ,x±_)t°- 



(2) 



with the generators of the fundamental representation of SU(3) and the 
classical color field of the nucleus in the color glass condensate description. 
In terms of the number density of color charges in the nucleus, pa{x_i_), its 
expression in the covariant gauge is : 

A+{x-,xi_) = -g6{x-)-^pa{x±) . (3) 

Turning this photon production amplitude into the amplitude for the pro- 
duction of a lepton pair, where the lepton has momentum ki and the anti-lepton 
has momentum k2 {ki + k2 = k), is straightforward. One simply needs to mul- 
tiply it by the free propagator of a virtual photon, which in the Feynman gauge 
is Dn^{k) = —iQuv/k'^, and by the leptonic current J'^ = ev{k2)^''u{ki). There- 
fore, the Drell-Yan amplitude is given by : 

M^^{p\q,ki,k2) = 2TTSip-- q-^ k-)Jd^x^e''-'^^+''^-P^>''^ [vix±)-l 

(4) 

It will be convenient to introduce a notation for the part of this amplitude 
that remains after we factor out the factor 2ttS{p~— q~— k~) : 

Mj,y{p\q,ki,k2) =2nS{p~- q~- k~)Mj,^{p,q;ki,k2) • (5) 

Indeed, as was shown in [34], the cross-section for the process qA — > ql~^l~A is 
then given in terms of M^^ by the following formula : 

1 d^qj_dq^ d^ki^dk^ d^k2±dk2 
'^"'^ " 2p- {2n)^2q- (27r)32fcf (27r)32fc^ 

x2Trdip- - q- - k^ - k^) ^ J2 (|M^,(p|q,fci,fc2)|')^ ,(6) 

*^ colors 
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where Eq, El, and E2 are respectively the energy of the final quark, of the 
lepton and of the anti-lepton. The bracket ( • • • )^ denotes the average over the 
color sources in the target nucleus. We have summed over the colors of the final 
quark and averaged over the colors of the incoming quark. When we sqiiare the 
amplitude in order to obtain the cross-section, it is customary to write it in the 
following way : 



(27r)4 M2 k{ 



(fki^dk-i d'^k2^dk2 _ ^2y^ j^^^ 



(7) 



where M is the invariant mass of the lepton pair (M^ = (fci + ^2)^), and Li^" 
is a leptonic tensor coming from the trace over the Dirac indices in the lepton 
loop^, 



1^1 1^2 



1^1 <^2 



ki ■ k2 

and where contains the factors that depend on the nuclear target : 



(8) 



X — tr 



(27r)2g 

V{x^)-l\ \v\v^_)-l 



S(p- -q- - k-) / d^x^d^y^ ^i{q^+k^)-(x^-y^) 



tr < 1^ 
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ffc)2 {p-kf 



(p-kf 



{q + kf 



(9) 



In this formula, we denote fc = fci + ^2, and we have assumed that the incoming 
quark has no transverse momentum: = 0. Of course, the integration over 
the variable q~ is trivially done thanks to the delta function. Note that the first 
trace is over color indices while the second one is over the indices carried by the 
Dirac matrices. 

By assuming translation invariance in the transverse plane for the nucleus, 
this formula may be rewritten in terms of the Fourier transform of the correlator 
of two Wilson lines : 



2p-{p 
xtr 



{q + kf {p-kf 

, 7*^(^ + 1^)7" 



ip-kf 



[q + kf 



where we denote 



^We also neglect the lepton mass. 



(10) 



nRlCil^) ^ I d^x^d'y^ e''-(---«-) i-tr (V {x ^)V\y ^)) ^ , (11) 
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being the radius of the nucleus. The function C{1±) is directly related to 
the Fourier transform of the cross-section for a qq dipole scattering on a target 
nucleus. Therefore, it could in principle be obtained from fits of Deep Inelastic 
Scattering on the same target. Note also that we have dropped terms with no 
Wilson lines, as well as terms with only one Wilson line. These terms have a 
6{q^+k±) under the integral, and correspond to a situation where no transverse 
momentum is exchanged between the target nucleus and the quark line. They 
do not contribute to the photon production cross-section. 

So far, we have written the Drell-Yan cross-section and structure functions 
for an incoming quark colliding on a nucleus. Naturally, they will be mea- 
sured for incoming nucleons. The relation between the quantities defined in 
eqs. (7) and (9) for the quark-nucleus system and the equivalent quantities for 
the proton-nucleus system are : 



dal^ = / dxi q{xi,ii^) da^^ 
Jo 



'p~ =Xiy/ s/2 



W^^ = f dx, q{xuH^) el W^\_^^^r^^ , (12) 

V 



where x\ is the longitudinal momentum fraction of the quark in the incoming 
nucleon. Note that the only combination of parton distributions of the proton 
that matter in this calculation is the structure function : 

F2(ari,M^) = 5Ze^g(a;i,M^), (13) 



since the "probe" with which we explore the proton is a photon. Most of our 
discussion will be about the structure functions of the quark-nucleus subsys- 
tem, and we will convolute them with the quark distribution of the proton (or 
deuteron in the case of RHIC) only at the end when we present numerical results. 



2.3 Structure functions and the Lam- Tung relation 

It is customary to perform a decomposition of the hadronic tensor W^^ into 
four structure functions. There are a priori three 4-momenta upon which W^^ 
might depend: P^^PJ^ the momenta of the proton and of a nucleon in the 
nucleus respectively, and k'^ the momentum of the virtual photon. This tensorial 
decomposition is usually performed in the rest frame of the lepton pair, where 
fc = 0. Denoting Ei and E2 the energies of the proton and of a nucleon in the 
nucleus in the rest frame of the lepton pair^, one first defines the following four 

^The vectors P^j should be kept in the center of momentum frame of the proton- nucleon 
subsystem, since it is in this frame that eq. (9) for W/^i, is defined. 
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projectors^ 



pnu 

TL 



pnv 



piiv 

L12 



P^P^ + P^P^ 



E1E2 



(14) 



^From these projectors, one defines the following four structure functions ; 



-'-TL — ^TL f^l^ ' 

y = pM^ w 



(15) 



Alternatively, different sets of structure functions can be defined by applying a 
4x4 "rotation matrix" to the previous set of structure functions. The helicity 
structure functions can be obtained in this way [41]. Several definitions of the 
helicity structure functions exist in the literature, and in this paper we are going 
to consider the so-called CoUins-Soper structure functions, defined as : 



T 



where the matrix R^^ is defined as : 



(16) 



Res - 2 



2 cos^ (5 



2 cos2 i 



- 



sin 25 



1 1+cos 5 
■ sin^ 25 



sin 25 



1+cos^ 5 
sin^ 25 



' 2 cos2 5 



1—3 cos 5 
sin^ 25 J 



(17) 



In this formula, we use the so-called "Collins-Soper frame" in order to define 

the angle 5: in the Icpton pair rest frame, the momenta of the two incoming 
hadrons are not opposite, and the Z-axis is chosen to bisect the angle between 
Pi and — P2. 5 \s the angle between the Z-axis and the momentum Pi, defined 
in figure (2) and given by 



cos 



2<5 



M2 



M2 + 



sm 



^5- 



M2 + fc 



(18) 



^Note that, since these projectors are defined in terms of the ratios / E\ and /E2, 
they are identical for the proton-nucleus system and for the quark-nucleus subsystem. This 
will become more obvious with the explicit formulas in eqs. (24), where one can see that 
the center of mass collision energy has disappeared. This means that we can perform these 
projections on the Wixu of the qA subsystem, and convolute with the quark distribution of 
the proton only later. 
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Figure 2: Definition of the angle 6. 



For the sake of completeness, we also show the expressions for the structure 

functions explicitly, 



W = -T 



2' 

2 , 1.2 



W = IM1±A[T -T ] 
w = _ .i^ 1 — ±11 1 — ±1 \T +T ] 

^^TL SM^kl 

(M^ + fci)(-2M^ + fci) 
At leading twist, the following relation 



= 0, (20) 

known as the Lam-Tung relation [39,40], is valid up to Next-to-Leading order, 
i.e. up to the order 0{as). 

2.4 Kinematics and matrix elements 

Below we will discuss violations of this relation in the kinematics appropriate 
to the RHIC and LHC experiments. Before that, it is instructive to be more 
specific. The vectors and are proportional to unit vectors along the light 
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cones* : 



pr = \r^vt, p,'' = ./j<, (21) 



2 

with 

""+ = 7!^^'°'°'^^' ^- = ^(^M-l). (22) 

(Remember that the proton is moving in the negative z direction and the nucleus 

in the positive z direction). The energies Ei and E2 of the proton and of a 
nucleon of the nucleus are, in the rest frame of the lepton pair, 

E^ = ^^s{M^ + kl), E2 = Q^s{M^ + kl), (23) 
so that the projectors of eqs. (14) can be rewritten as : 

Performing exphcitly the Dirac algebra and the contractions with these pro- 
jectors, we obtain the following expression for the four T's : 

where a e {TL, LI, L2, L12} with 

T -iJ (l-^)^i I [zQ±-{l-z)k^? 
[ [M2 - 2p • ^ (1 - ^)[M2 + 2g • 

[(1 - z)gi + (2 - z)g^ ■ fc J 1 
[M"^ - 2p ■ k][M'^ + 2q ■ k] J' 

^ 32(1-0) , . r 1 1-z 



M'^ + 2q-k M'^-2p-k 
32z2 M2 r9i(gi + fei)2 



"^^^ 1-z (M2 + fel)2\ [M2 + 2(7-fc]2 

T -61 g^-(g^ + fc±) [ 1-z 1 

M^ + k'i [M2 + 2g • fc] \ M2 - 2p • /c M2 + 2g • /e 

(26) 

*We can aiTord to be a bit sloppy here about the use of the squared center of mass energy 
for the proton- nucleus system, s, instead of that for the quark-nucleus subsystem, s = xis. 
Indeed, since the energies Ei and E2 are also proportional to ^/s, this choice is irrelevant in 
the definition of the projectors. 
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Eqs. (25) and (26) thus provide the full leading result for the production of a 
lepton pair in proton-nucleus collisions in the CGC framework. These results 
are for a quark with a fixed longitudinal momentum p~ . In order to obtain the 
corresponding results for an incoming proton, one needs simply to convolute 
with the quark distribution q{xi), and use eqs. (26) with z = k~ /{xi\/ s/2). 

This result is of the same order in as the standard NLO qg — > 57* — > 
q£^£~ contribution, but compared to the standard pQCD result it resums all 
the multiple scatterings on the dense nucleus. One point should be emphasized: 
our calculation incorporates also the "LO" process (79 — > 7* — > £~^l~ , although 
in a non obvious way. Indeed, the intermediate virtual and space-like quark 
in the diagrams of figure (1) must be interpreted in some kinematical domain 
as an antiquark belonging to the wave-function of the nucleus (in that case, 
the quark in the final state came also from the nucleus wave- function) . Note 
however the pQCD NLO process qq — » gj* — » gi~^£~ would only appear at the 
next order in our calculation, via diagrams similar to those of figure (1) but with 
a gluon in the final state in addition to the photon. In order to understand why 
this process does not appear on the same footing here, one has to remember 
that in our approach the antiquarks are produced explicitly from the gluons, 
and that a process that has both an antiquark in the initial state and a gluon 
in the final state is non-leading in our framework. 

2.5 The collinear factorization limit 

In eq. (10), the function C{qj_ + k±) resums all the multiple scatterings that 
the quark undergoes while it goes through the nucleus. However, it is possible 
to extract from this formula the contribution from a single scattering, i.e. the 
leading- twist contribution. In this limit, one should recover the results from 
collinear factorization in perturbative QCD. 

This is done as follows. One assumes that the momentum scales qj_,k±, M 
that characterize the final state are much larger than the saturation scale, and 
one makes the approximation that the transverse momentum of the incoming 
gluons, l± = qj_ + k±, is large compared to Qa but small compared to the scales 
of the final state. This approximation is legitimate in this regime, because when 
there is a large range between Qg and the scales of the final state, the integration 
over Z_L is dominated by a large logarithmic contribution that comes from values 
of /_L comprised in this range. One must first approximate the matrix elements 
Ta by assuming that the transverse momentum of the gluons coming from the 
nucleus, l± = + k±, is smaller than the other scales. As pointed out after 
eq. (10), these matrix elements vanish when l± = 0, and therefore one should 
keep higher orders in the expansion in Formally, one can write : 

T„ = + (27) 

Since the function C{1±) depends only on the modulus of the vector /j_ but 

not on its orientation, the term linear in l± in this expansion will vanish when 
inserted into eq. (25) when we perform the integration over the orientation of 
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the vector l±. We therefore need to expand the to second order in l±. 
Assuming this expansion has been performed, we can now write 

RecalUng the fact that in the leading-twist hmit - vahd here since lj_ Qg -, 
l^C{l±) is proportional to the non- integrated gluon distribution of the nucleus, 
the integral over l±_ that appears in the previous equation gives the usual inte- 
grated gluon distribution, xG^ (x, Q^), with a resolution scale which is the 
upper limit of the range of validity of the expansion in eq. (27). 

It is straightforward to perform explicitly the expansion in powers of l±_ of 
the matrix elements 7^, in order to obtain the second order coefficients U^*. 
One gets : 



=64^^(1 -2) 



[ki + (1 - Z)M2]4 
M2 kl 



(ki + M2)2 [ki + (1 - 0)M2]2 



B[\,=-128^^(l-^) 



3n ^2 Af2fci 



' (fei + M2) [kl + (1 - 0)M2]3 ' 

=32^2(1-^)^^ 1 

^ [fci + (1 - Z)M2]4 

X [4(1 - zfklM^ + [1 + (1 - zf][k'i + (1 - zfM^W . (29) 

Using these expressions, it is easy to verify that the Lam- Tung sum rule given 
in eq. (20) holds identically in the leading twist (single scattering) limit. 

2.6 Angulcir coefficients 

Based on the general Lorentz structure of the cross section, one can write the 

transverse momentum integrated Drell-Yan cross section in terms of the polar 
and azimuthal angles of the dilepton pair in the dilepton center of mass frame 
as 



dN 

'dn 



da 



da 



dm'^k 



(30) 



where dil = dcosOdfj) is the solid angle, with 9 and </> defined in Fig. (2) and k 
is the dilepton four momentum with 



dN 



W^^ (1 + cos^ 6*) -h ^M^^ (1 - 3 cos^ 6*) 



sin 2e cos (j) + W^^ siv? 6 cos 2^ 



(31) 
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It is customary to parameterize this cross-section in terms of the angular coef- 
ficients Ao,Ai,A2 as follows, 



dN 



3 
16^ 



1 + cos^ e + -Ao (1-3 cos^ e) 

+Ai sin 26 cos + sin^ 9 cos 2(j) 



so that 



^0 



Ai 



2W. 



(32) 



(33) 



in terms of which the Lam- Tung relation simply reads Aq = A2. Below we 
will show our results for the structure functions W's as well as the angular 
parameters A's. 



3 The dipole cross-section 
3.1 Balitsky-Kovchegov equation 

All single inclusive particle production cross sections in the Color Glass Con- 
densate formalism depend on the Fourier transform C{1±) of the dipole cross 
section^, which is the probability for scattering of a high energy quark-antiquark 
dipole on a target. The dipole cross section evolves with rapidity according to 
the JIMWLK equation, so that given the dipole profile at some initial point 
xo, the dependence of the cross section on x (or rapidity) can be in principle 
obtained from the JIMWLK equation. In practice, this is very difficult to do 
(see [42]) and it is much simpler to use an approximate form of this evolution - 
known as the BK equation [29,30] - which can be obtained from the JIMWLK 
equation as a mean field approximation in the large limit. Here we will use 
the numerical solution of the BK equation in order to investigate the properties 
of the Drell-Yan structure functions. 

In momentum space, the BK equation takes the form 

dT{kpy) =as{K(^T){k^,Y)-asT\k^,Y) , (34) 

where we denote ctg = agNc/n. The operator K is the well known BFKL kernel 
in momentum space [43,44], whose action on the function T is given by 

+00 ^ 

dik'l) J k'lTik'^,Y) - klT{k^,Y) ^ klTik^, Y) 



2 ( 

{K^T){k^,Y)^ I ^ 



H-k'll ^4k'l + kl 



(35) 



^Note that in case of gluon scattering one has a dipole in the adjoint representation. 
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The function T{k±, Y) is the Bessel- Fourier transform of the dipole-target scat- 
tering ampUtude T[r±, Y): 

+ CXD 

T{k^,Y)= j ^Joik±r±)Tir±,Y) , (36) 



where r± is the size of the qq dipole and k± is its conjugate transverse momen- 
tum. The dipole ampHtude T is defined in terms of the correlator of two Wilson 
lines of gauge fields in the target as 

T{ru y) = 1 - i-TV (yt(O)y(r^))^ , (37) 

where we have assumed translation invariance in the transverse plane in order to 
set the quark transverse coordinate to 0. The scattering amplitude T is related 
to the function C{1±) by 



C{1^) = Jcfr^e''^-^^ [l-T{r^,Y)] 



(38) 



Therefore, in order to obtain the function C{1±) needed in the evaluation of 
the Drell-Yan structure functions, we must first solve the BK equation 34 for 
T{k±,Y), then invert eq. (36) in order to obtain T{r±,Y), and finally obtain 
C{1±) by a Fourier transform. 

The BK equation must be siipplementcd by an initial condition at 1" = 0, for 
which we take the value of T{k±) predicted by the MV model [45-47] with an 
initial saturation momentum of = 2 GeV^. This initial condition is assumed 
to correspond to an X2 in the nucleus of X2 = 10~^, which means that the 
parameter Y in T{r±, Y) should be interpreted as F = ln(a;2/a;2). For values of 
X2 larger than a;^, we use the following naive extrapolation : 

C(i±,a;2)= C(/i,4). (39) 

The idea behind this extrapolation is that the unintegrated gluon distribution 

that appears in the gluon yield (see [18]) is proportional to C{1±_), and should 
vanish like (1 — a;2)'* when X2 approaches 1 (the exponent 4 comes from quark 
counting rules). 

Note also that we solve the BK equation with a fixed value of the coupling 
constant a^. However, we set it to a rather low value in order to correctly 
reproduce the rate of growth of the saturation scale, as inferred for instance 
from the study of geometrical scaling at HERA [48] . This approach has already 
led to good results in the analysis of limiting fragmentation in [49] . 

3.2 Models of the dipole cross section 

As an alternative to solving the BK equation, one can model the behavior of 
the dipole cross section based on its known properties. These phenomenological 
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models have been used to investigate particle production at RHIC and LHC 
as well as structure functions in Deeply Inelastic Scattering in HERA. Here we 
will use two of these models, known as the KKT [7] and DHJ [9,10] parame- 
terizations, to calculate rapidity and transverse momentum dependence of the 
Drell-Yan structure functions in the kinematics appropriate for RHIC and LHC. 

The KKT parameterization of the dipole cross section was quite successful in 
description of the forward rapidity hadron production cross section in deuteron- 
gold collisions at RHIC. In this parameterization, the dipole forward amplitude 
is modeled as 



T{r±,Y) = 1-exp 



(40) 



In the KKT parameterization, the anomalous dimension 'y{r±,Y) is given by 

,(.., y) = 1 f 1 + ] (41) 

and 

^' ' (A/2)(F-Fo) ^ ' 

The saturation scale is given by (5,5(y) = Qo cxp[A(y — yo)/2] with y = In l/a:;2 
and Yq = 0.6, A = 0.3. While the KKT parameterization of the dipole cross 
section was the first one to reproduce quantitatively the forward rapidity RHIC 
data, it has some shortcomings. For example, the assumed dependence of the 
anomalous dimension on r±_ (or equivalently k±) is too flat as shown in [9]. 

In the DHJ parameterization [9,10] of the dipole cross section, the anomalous 
dimension is given instead by 

7(r^, y) ^ 7. + (1 - 7.) ^^liiZl^I^ (43) 

' ^''xY + logil/rlQ^,)+dVY ^ ' 

with 7s ~ 0.62, and d = 1.2. This model has the advantage that it accounts 
for the more rapid change of the anomalous dimension with transverse size 
(or momentum) which is needed in order to quantitatively describe the hadron 
transverse momentum spectra in both forward and mid rapidity deutcron-gold 
collisions at RHIC. For further details of the two models and their applications 
to RHIC data, we refer the reader to [9,10]. In the following we do not consider 
the KKT model any further, and we use either the BK evolution equation or 
the DHJ model. 



14 



4 Numerical results 



4.1 Kinematics 

Let us complete the eqs. (26) by some useful formulas for the denominators that 
appear in the 7^'s : 

z 

,.2 ^ , M"^ {zq,-il-z)kA_f 

^ +^^-^ = -+ .(1-.) ' ^''^ 

where we recall that the variable z is defined as 2; = k~ /p~ . These formulas 
highlight the fact that the 7^'s defined in eqs. (26) depend on the parameters 
M, fcj^, y of the virtual photon, on the transverse momentum of the outgoing 
quark, and on the momentum fraction z of the photon relative to the incoming 
quark. In particular, this variable z is the only reference to the momentum of 
the incoming quark in these formulas. 

For an incoming quark carrying the fraction a;i of the longitudinal momen- 
tum of the proton, we also have : 



zxi = d — - e y . (45) 

Note that since both z and xi must be smaller than unity, this implies obvious 
limits to the phase-space allowed for the produced virtual photon, since we must 
have : 

< 1 . (46) 

s 

The notation for this quantity comes from the fact that it is also the minimal 
allowed value of xi. Thus, the structure function Ta for the proton- nucleus 
system can be written as 

1 

Tt= J dx, g(x„M^) "^l^ J ^Ciq,+k^,x,)T„(z 

(47) 

One last point remains to be clarified: the dependence of the function C on X2 
is now explicitly shown. This dependence arises from the renormalization group 
resummation of small x effects via JIMWLK equations and has the interpreta- 
tion as the longitudinal momentum fraction at which the nucleus is probed. It 
can be evaluated as 

X2 = ^^^^ = <"e'^ + , . (48) 

Note that the minimal value of xi is 



X\ 



^min = W + , (49) 
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(With our choice to have the nucleus moving in the +z direction, small values 
of X2 are reached when the photon rapidity is negative.) and that the relation 
between xi and X2 is more symmetric if written as : 

{xi - xf'^Xxz - ") = ^ . (50) 

Moreover, we can if necessary replace the integration over xi in eq. (47) by an 
integration over X2 by noting that 

dTPA _a;i-a;f" dT^^ ^^^^ 



dx2<Pq±_ X2 — a;™ dxi(Pq^_ 

This is useful in order to determine what are the values of X2 that contribute 
the most to the Drell-Yan structure functions. 



4.2 xi and X2 dependence of the integrand 

Let us first study the X\ and X2 dependence of the integrand. We do this for 
pA collisions at the LHC center of mass energy, i.e. = 8.8 TcV, for nuclei 
with A^^^ = 6. In figure 3, we display the Xi dependence for the four W's (i.e. 
we have multiplied eq. (47) by the rotation matrix i?^, , given in eq. (17)), for 
fc_L = M = 3 GeV and y = -3. We use the CTEQ6 (sec [50] for details) set of 
parton distributions for the quark and antiquark distributions of the proton. In 
this figure, we compare the x\ dependence that results from the BK evolution 
of the dipole cross-section, and from the DHJ model. As one can see, the two 
models lead to results that are very similar: the largest contributions come 
from values of x\ immediately above the minimal value a;™'", and the integrand 
drops very quickly as xi increases. One can also note that the integrand for 
vanishes and changes sign at some intermediate value of xi (for this integrand, 
we have displayed its absolute value, hence the cusp at the location of the zero 
on the logarithmic plot). In fact, for the kinematical parameters of the virtual 
photon used here, the two descriptions of the dipole cross-sections are very 
close even quantitatively. However, as we shall see later, this is not the case 
everywhere. 

In figure 4, we have used the eq. (51) in order to display the X2 dependence 
of the same integrands. In order to see how the dominant values of X2 vary with 
the photon rapidity, we have considered three values oi y, y = —6, —3 and 
(the first one can be reached only by the CMS and ATLAS detectors, while the 
last two may also be studied by ALICE). To avoid overcrowding the figure, we 
display the results of this study only in the case of the BK evolution^. As in 
the case of the Xi dependence, the dominant values of X2 are located just above 
the minimal value a;™™, and the integrand decreases quickly as X2 increases. 

^At X2 = 10^^, one can notice a very small discontinuity in the slope of all the curves. 
This is the value of X2 where we switch to the extrapolation described by eq. (39). As long 
as this point is far from x™'", this extrapolation will not affect significantly the result of the 
integration over X2- 
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Figure 3: xi dependence of the integrand for the W's. The kinematical param- 
eters of the virtual photons are fcj_ = 3 GeV, M = 3 GeV and y = —3, and 
the center of mass energy is y/s = 8.8 TeV. Lines: BK evolution. Points: DHJ 
model. 



This enables one to select different windows in X2 by varying the kinematical 
parameters of the virtual photon. This is illustrated in figure 4 with three 
different values of the photon rapidity. As one can see, it is possible to explore 
very small values of X2 by going to negative rapidities, provided that the photon 
transverse momentum and invariant mass remain rather small. 

4.3 Structure functions for pA and dA collisions 
4.3.1 k± dependence 

Let us now display the results for the W's, integrated over the momentum 
fraction xi of the incoming quark. We have performed this calculation both for 
the LHC and for RHIC. For the latter, we have used a center of mass energy 
y/s = 200 GeV and the incoming projectile is taken to be a deuteron^ rather 

'^The CTEQ package does not provide parton distributions for the deuteron. In order 
to overcome this restriction, we assume that one can just add the parton distributions of a 
proton and of a neutron, i.e. we neglect the shadowing that may occur in the deuteron. This 
is a rather good assumption as long as xi is not too small. In order to obtain the quark 
distributions in the neutron, we just exchange the u and d valence distributions of a proton, 
and we assume that the sea quark distributions are identical in a proton and a neutron. 
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LHC - k = 3 GeV , M = 3 GeV , Y = -6,-3,0 
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Figure 4: X2 dependence of the integrand for the W's, with the BK evolution. 
The kinematical parameters of the virtual photons are A;_l = 3 GeV, M = 3 GeV 
and we consider the rapidities y — —6, —3 and 0. The center of mass energy is 
^/i = 8.8 TeV. 



than a proton. 

In figure 5, we display the dependence of the Drell-Yan structure functions 
on the transverse momentum of the virtual photon, for a fixed invariant mass 
and rapidity (this rapidity is chosen negative to be in the region where the 
nucleus is explored at small X2). The calculation has been performed with the 
function C{l^,X2) resulting from the BK evolution (lines) and from the DHJ 
model (points), both for the LHC (left) and for RHIC (right). We also display 
on the same plot the difference W^^ — which is expected to vanish if 

the Lam- Tung relation is satisfied. Since we have checked explicitly that this 
combination vanishes identically in the leading twist approximation, any non- 
zero value for this quantity is due to higher twist corrections ~ i.e. to multiple 
scatterings of the quark in the nuclear target. One can see in figure 5 that the BK 
evolution and the DHJ model lead to qualitatively similar features in the Drell- 
Yan structure functions. The difference between the two descriptions is never 
very large (at most a factor 3), but one can see that it is more pronounced at 
large k±, and also more pronounced at RHIC energy compared to LHC energy. 
Moreover, the fact that the W^'s obtained with the DHJ model drop slightly 
faster at large k± than those obtained from the BK evolution suggests that the 
difference comes from the anomalous dimension 7. In the BK evolution, it is 
always slightly above 1/2, while the DHJ parameterization is such that this 
anomalous dimension goes back to unity at large momentum. 

Note that the difference W^^ — should be compared to the structure 
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Figure 5: Drell-Yan structure functions as a function of k±, for RHIC (right) 
and LHC (left). Lines: BK evolution. Points: DHJ model. 



function W^^^ (indeed, the W^^,^ structure function, obtained by contracting with 
gf^", is closely related to the total cross-section). One sees that this combination 
is more and more suppressed as the transverse momentum of the virtual photon 
increases, as expected for a higher twist quantity. This means that violations of 
the Lam- Tung relation should be searched at low k± rather than at large k± . In 
figure 6, we display the angular coefficients defined in section 2.6 for the same 
set of parameters. Violations of the Lam- Tung relation should are visible in the 
fact that the difference Aq — Ai is non-zero. One can see that they decrease, 
albeit very slowly, when the transverse momentum increases. 

4.3.2 M dependence 

We display in figures 7 and 8 the same structure functions and angular coeffi- 
cients as a function of the virtual photon invariant mass. Again, the calculation 
was done both for RHIC and the LHC, with the BK evolution and the DHJ 
model. We observe for the M dependence results that are qualitatively similar 
to what we already found for the fc^ dependence. The DHJ model and the BK 
evolution are closer at the LHC rather than at RHIC. And again, we see that 
the suppression of the difference W — W ^ /2 compared to the "total" struc- 
ture function W^^ is larger at large mass, or conversely that violations of the 
Lam- Tung relation are larger at small invariant mass. 

4.3.3 Y dependence 

Finally, we present results concerning the rapidity dependence, at fixed and 
M, of the Drell-Yan structure functions and angular coefficients in figures 9 and 
10. Note that the vertical scale in figure 9 is now linear, in order to better show 
the variations with Y and the differences between the BK evolution and the 
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Figure 6: Angular coefficients as a function of k±, for RHIC (right) and LHC 
(left). Lines: BK evolution. Points: DHJ model. 



DHJ model. This natural makes the differences between the DHJ model and 
the BK evolution look larger than in the previous plots, while they are in fact 
comparable in magnitude. 

Now, one sees qualitative differences between the two descriptions of the 
nuclear target. In absolute terms, the two descriptions do not differ by vast 
amounts, but the variations with Y are quite different. The general trend is 
that the BK evolution (with the parameters described in section 3.1) leads to 
significantly stronger variations with the rapidity. One should keep in mind 
that the DHJ model has been primarily devised to reproduce correctly the 
k± dependence of the hadronic spectra observed by the BRAHMS and STAR 
collaborations at RHIC, and in particular their slope and how this slope varies 
with rapidity. However, it did not predict the absolute normalization of these 
spectra, and a "K-factor" had to be adjusted by hand. Moreover, the K-factor 
necessary to reproduce the normalization of the data turned out to be different 
at each rapidity. What we observe here may be related to this fact, and seems to 
indicate that the DHJ model should still be improved in order to get correctly 
the changes in the absolute normalization with rapidity. 

Regardless of the description used for the function C{l±,X2) - BK evolution 
or DHJ model -, one also sees very different shapes at RHIC and the LHC. In 
the DHJ model, the slope tends to be inverted from RHIC to the LHC. And 
with the BK evolution, one goes from fairly symmetric shapes at RHIC to a 
very asymmetric rapidity distribution at the LHC. However, one should take this 
remark with caution, because our calculation of the Drell-Yan structure function 
treats the proton and the nucleus in very different ways, and is only valid when 
X2 is small in the nucleus and xi is rather large in the proton. Therefore, the 
accuracy of our results at positive rapidity may suffer significantly from this 
treatment. 

Finally, one can note that the violations of the Lam- Tung relation, that one 
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can see in figure 10 on the difference Aq — A2, rise steadily when the photon 
rapidity goes from positive to negative values (at least if one excludes the frag- 
mentation regions of the two projectiles). This is also expected of a higher twist 
effect: indeed, at smaller y, i.e. at smaller X2, the density of color charges in 
the nucleus increases, which leads to more important rescattering effects. 



5 Conclusions 

In this paper, we have studied in the Color Glass Condensate framework the 
production of Drell-Yan pairs in proton-nucleus and deuteron-nucleus collisions 
at RHIC and at the LHC. Our treatment is asymmetric, because we only use 
the CGC description for the nucleus, while the proton or deuteron is described 
by means of the usual parton distributions. 

All the nuclear effects that affect the Drell-Yan structure functions are en- 
coded in the function C{1±), which is the Fourier transform of the dipole- nucleus 
cross-section, and we have considered two descriptions of this object: the DHJ 
model and the solution of the BK equation (with the McLerran-Venugopalan 
model as the initial condition). Overall, the most striking differences between 
these two description are seen in the rapidity dependence of the structure func- 
tions, where the BK description leads to a much stronger rapidity dependence. 

We have also discussed the Lam- Tung relation, a linear combination of the 
structure functions which is known to vanish identically up to (and including) 
the Next-to-Leading order in the leading twist approximation. Such a quantity 
is indeed very interesting to investigate in order to study the effects of higher 
twist contributions, because any non-zero value is by definition coming from 
higher twist. We observed for this quantity a behavior which is indeed what one 
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Figure 8: Angular coefficients as a function of M, for RHIC (right) and LHC 
(left). Lines: BK evolution. Points: DHJ model. 

expects from a higher twist quantity: it drops at high momentum or high mass 
slightly faster than the total cross-section, and it is larger in regions of rapidity 
where the nucleus is probed at smaller a;, i.e. where it has a larger parton 
density. A detailed experimental study of Lam- Tung relation then should prove 
very fruitful in advancing our knowledge of QCD in its novel high parton density 
regime. 
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